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Abstract 

Until recently, the preferred strategy to identify stationary points in the scalar 
potential of S f O(8)-gauged Kf — 8 supergravity in D — 4 has been to consider 
truncations of the potential to sub- manifolds of E 7 ( +7 )/SU(8) that are in- 
variant under some postulated residual gauge group G C £0(8). As powerful 
OXJ' alternative strategies have been shown to exist that allow one to go far beyond 

what this method can achieve - and in particular have produced numerous 
solutions that break the SO (8) gauge group to no continuous residual sym- 
metry - independent verification of results becomes a problem due to both 
the complexity of the scalar potential and the large number of new solutions. 
This article introduces a conceptually simple self-contained piece of computer 
code that allows independent numerical validation of claims on the locations 
of newly discovered stationary points. 
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£h! 1 Introduction 

& . 

Maximally extended (Af = 8) supergravity in four dimensions recently has seen a 
resurgence of interest, in particular due to speculations about its potential per- 
C*~) ■ turbative finiteness [Hll||, as well as the AdS4/CFT3 correspondence in the 

context of ABJM/BLG theory [4] and also early ideas about an AdS/CMT cor- 
respondence mis]. 

As is well known, Af = 8 supergravity allows the promotion of its 28 vector 
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fields to nonabelian gauge fields. This was first demonstrated for the gauge group 
SO (8) in [J] . As has been shown in [JJJ[TU], there are other, more exotic, options as 
f- ^ ■ well; in particular, it is possible here to introduce noncompact gauge groups such 

as S0(p,8-p). 

In order to maintain supersymmetry, such a deformation of the model requires 
the introduction of a potential on the 70-dimcnsional manifold of scalars. This 
potential is known to have a number of nontrivial stationary points, some of which 
correspond to stable vacuum states with broken symmetry and usually some residual 
jrt ■ supersymmetry. Unfortunately, the analysis of these potentials is fairly involved, as 

they are most naturally regarded as functions on the coset manifold E 7+7 ^/SU(8), 
which is infcasiblc to parametrize analytically by introducing coordinates, as is 
readily seen by simple order-of-magnitude guesstimations on the number of expected 
combinations of trigonometric and hyperbolic function factors in a generalized Eulcr 
angle parametrization. 

For a long time, the favored approach to nevertheless extract some information 
about stationary points was to consider restrictions of the potential to sub-mainfolds 
of E 7+7 - ) /SU(8) that are invariant under some subgroup G of the gauge group, the 
rationale being that the stationarity condition on the restricted manifold readily is 
lifted to the full manifold, as the first order term in a perturbative expansion of the 
potential must be a G-invariant scalar. While further technical subtleties such as 
coordinate singularities might complicate the analysis, this has been shown to be a 
powerful technique that allowed the determination of all S'C/(3)-invariant stationary 
points [TT], as well as a further stationary point with residual symmetry £0(3) x 



50(3), cf. [12]. No further solutions have been found until a new strategy has been 
proposed in |13j that increased the number of known nontrivial solutions from six 
to twenty. 

While this method, which is based on numerical techniques that were discovered 
at about the same time [14] as the scalar potentials of gauged maximal supergrav- 
ity, allows a much deeper analysis and is expected to give rise to numerous other 
solutions - as it did [15] for 5*0(8) x 5*0(8)-gaugcd Chcrn-Simons Af — 16 super- 
gravity in D = 3 |16j - it raises two new issues: First, this is 'only' a numerical 
approach that provides evidence for new stationary points, but no stringent proof of 
their existence. And second, due to the large number of solutions and the analytic 
complexity of these potentials, checking these claims unfortunately amounts to a 
substantial task. 

The first issue can be addressed in a semi-automatic way by employing another 
modern and fairly revolutionary algorithm - the the PSLQ 'lattice shift' algorithm 
to find integer relations between numbers given with high-precision |17j . The basic 
principle is demonstrated in section 7 of |18j and also briefly explained in section [3] 
As obtaining a very large number of valid digits is computationally expensive, 
the cut-off parameters for the PSLQ algorithm will generally require some man- 
ual tweaking, so striving for full automatization of this step, while possible, would 
probably not be a worthwhile objective. The post-processing of analytic conjectures 
obtained via the PSLQ algorithm, which amounts to establishing their analytic va- 
lidity, should be fairly straightforward to automatize, however. The challenge is to 
analytically exponentiate exact expressions for Lie algebra elements using symbolic 
algebra (probably somewhat tailored to the algebraic task) and then symbolically 
verifying the stationarity condition. 

The second issue is addressed by this article, which provides computer code to 
numerically validate claims about new stationary points of S r O(8)-gauged Af = 8 
supergravity with little effort. As it is expected that many more solutions and 
maybe also other, even more powerful methods to analyze such scalar potentials 
could be discovered in the future, having a self-contained dedicated tool for such a 
task - as is presented in this article - seems to make sense. 

2 The Code 

Given the main objective of validation of stationary points in the scalar potential 
of SO (8)-gauged Af = 8 supergravity, computer code that calculates this potential 
should be clean, compact, self-contained, and easily verifiable. This is best achieved 
by an independent re-implementation of numerical code that emphasizes simplicity 
and verifiability (i.e. a close match between mathematical formulae and code) 
over algorithmic tricks to achieve high performance and also is independent of the 
codebase that was used to find the new solutions presented in 



2.1 Design Choices 

The problem of numerically computing the supergravity scalar potential - while 
involving some intricate definitions - only requires arithmetics, a small bit of com- 
binatorics, and as the only tricky step, exponentiation of a complex matrix. There- 
fore, pretty much any programming language could be used as an implementation 
basis. While problems of such a nature that do not require symbolic transformations 
often are addressed with MATLAB [19], or its free alternative GNU Octave [20] . 
Python [21] has been chosen for this work, for the following reasons: 



1 That codebase is at the time of this writing a complex multi-language mix of LISP and Python 
modules that mostly deal with problems other than validating results and are not even available 
on all computing platforms (in particular Microsoft Windows). It will be provided after major 
clean-up work to make it easy to install also on non-Unix platforms. 



• Python allows a transcription of the required combinatorics to computer code 
with considerably less overhead than many other languages (but not as effec- 
tively as Lisp or Scheme). 

• Python is free (in contradistinction to commercial computer algebra packages) 
and readily available across a broad range of hardware and operating system 
platforms. 

• Python syntax is very simple and conceptually sufficiently similar to other 
programming languages popular with casual programmers that it can be rea- 
sonably expected to be comprehensible to a broader audience than many other 
languages. 

• Python provides a command prompt that allows interactive evaluation of code. 

• Python itself has become somewhat popular for many scientific problems. 

A major architectural drawback of the code presented here is that, as it de- 
liberately is kept algorithmically simple, it neither uses sophisticated combinatorics 
nor efficient black-box approaches (such as relational database algorithms for sparse 
tensors, as the LambdaTcnsor [32] package does) to keep the computational effort 
to a minimum. However, for the purpose of validating results, speed is not an 
important concern. 

2.2 Installation and Usage 

The codebase utilizes the Python modules numpy for (non-sparse) tensor numerics 
as well as scipy (Scientific Python) for matrix inversion and exponentiation, and 
hence requires both Python as well as these extension modules to be installed. Most 
free Unix distributions are strongly modularized and provide the corresponding 
packages; for other platforms such as Microsoft Windows, convenient pre-bundlcd 
Python packages that contain these modules are available. Installation of Interactive 
Python (ipython) also is highly recommended. 

Once Python is installed and the source that accompanies this article is down- 
loaded from http://arxiv.org/e-print/1007.0600 and unpacked, the functions 
described below can be made available by importing the corresponding module into 
Python. The most important constant certainly is the dictionary vacua_v70 from 
e7_vacua.py - which contains all the known stationary points. The dictionary keys 
are strings of the form '#0', '#1', '#2', etc. that correspond to the numbers used 
in the tables in [T3]. All Python functions come with interactive online help. A 
transcript of an interactive Python session is given in figure [1] 

For some investigations, scripting is preferable to interactive usage. This is most 
easily done by adding new Python files to the code directory; alternatively, one can 
extend the PYTH0NPATH environment variable to tell Python where to look for mod- 
ules and put additional code into an arbitrary directory. The file e7_example . py 
which also is listed in figure [5J shows scripting use. It is executed as 'python 
e7_example .py'. 

2.3 Modules and Functions 

The code itself is modularized and consists of these components: 

2.3.1 The tensor_io.py module 

This component defines functions that read and write numerical (higher-rank) ten- 
sors in a well-defined simple sparse textual data format that is also human-readable, 
other codebases. It is generally expected that programs which manipulate tensors 
might want to save these to files, both for persistent data storage and data exchange. 
The full format definition, which may be adopted by other programs, is given in the 
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ipython 








In 


[1] : from e7_analysis import 


* 






# 


(This takes a while, as definitions have 


to 


be computed) 


In 


[2] : from e7_vacua import * 








In 


[3]: help ( e7_Al_eigenvalues ) 








# 


This brings up online help for this fund 


ion 


In 


[4]: e7_Al_eigenvalues (vacua 


_v70 ["#3"]) 






Out [4] : 








(array([ 1 . 33758921 -0 . 09352799 j 


1.33758921 -0.09352799 j , 




1 .09478 134 + 0. 00705399 j , 


1.33758921 


-0. 


09352799 j , 




1.33758921 -0.09352799 j , 


1.33758921 


-0. 


09352799 j , 




1.33758921 -0.09352799 j , 


1.33758921 


-0. 


09352799 j]) , 


1 


!l. 4999987205717096 , 










1.4999987205717116 , 










1.0000000000469873 , 










1 .4999987205717116 , 










1.4999987205717089 , 










1.4999987205717109 , 










1 .4999987205717116 , 










1 .4999987205717105] ) 








IB 


[5]: e7_test_stat ionar ity_Q ( 


vacua_v70 ["SU(4) "]) 


Out [5] : 3.0086481280021352e-15 








# 


This shows that the Q- tensor 


stati onari ty 


condition is 


# 


satisfied to high numerical accuracy for 


the known 


# 


stationary point with SU(4) s 


ymmetry . 







Figure 1: An interactive ipython session with the validation code. 



Python documentation of this module. It provides the functions tensor_print(), 
tensor jreadO , and tensor_write() . 

2.3.2 The e7_combinatorics .py module 

This component provides a few combinatorics-related definitions that are of use to 
other components. In particular, it defines two-index and four-index ranges running 
from (0,0) to (7,7) and (0,0,0,0) to (7,7,7,7), respectively. Also, it introduces 
iteration over 3-permutations and 4-permutations. Here, it has to be kept in mind 
that index counting starts at 0, not 1. Matching published data to numerical results 
hence requires shifting all index ranges found in the literature by one. 

2.3.3 The e7_def initions .py module 

This component provides definitions related to the spin(8) algebra that are then 
used to define explicit forms for the 133 complex 56 x 56 E 7 r +7 \ generator matrices 
of the fundamental representation. The most important definitions are these non- 
sparse numpy higher-rank arrays: 

• T_e7 - the 133 x 56 x 56 tensor T^ E7 \/ from formula (A12) in [TB] , 

• so8_gamma_a_aS_aC - the 8x8x8 tensor r y a a a containing spin(8) 'gamma 
matrices'. The order of indices is vector, spinor, co-spinor, as indicated by 
the name. 

• so8_gamma_ab_aS_bS, so8_gamma_ab_aC_bC, 
so8_gamma_abcd_aS_bS, so8_gamma_abcd_aS_bS - likewise, these (non-sparse) 
8x8x8x8 and 8x8x8x8x8x8 arrays give the tensors 7 a6 Q ^, 7 at ' & g, 

„.abcd nT1 j „,abcd 

7 a 0, and 7 A/ j. 



from e7_vacua import * 
from e7_analysis import * 

for n in range (21): 
key= "#"+str (n) 
v70 = vacua_v70 [key] 
(gens , svd)=e7_residual_gauge_group_generators (v70) 

print "Solution:", key , \ 

"P=" ,e7_potential_from_v70(v70) ,\ 
"Q=" ,e7_test_stationarity_Q(v70) ,\ 
"dim(res. GG) = " , len (gens) 

-6.0 Q= 0.0 dim(res. GG) = 28 

-6.68740304976 Q= 2. 12995452535e-15 dim(res. 

-6. 98771242969 Q= 5 . 07925733935e -15 dim(res. 

-7. 19157562605 Q= 6.61364975244^-06 dim(res. 

-7.79422863406 Q= 3.98414774796&-15 dim(res. 

-8.0 0= 3.008648128e-15 dim(res. GG)= 15 

-14.0 Q = 2.98299102544&-15 dim(res. GG)= 6 

-9.98708340034 Q= 1 . 69023258423e-06 dim(res. 

-10.434712595 Q= 4. 19318631751e-08 dim(res. 

-10 .6747542478 Q= 7.60840475615e-05 dim(res. 
-11 . 6568542495 Q= 3. 12755317279e-07 dim(res. 
-11.9999999998 Q= 8 . 08047489882e -07 dim(res. 
-13.6236525917 Q= 1 . 258980691 le -06 dim(res. 
-13.6761142184 Q= 5.80889427017^-07 dim(res. 
-14.97038467 Q= 1 . 4218350389e -07 dim(res. GG)= 1 
-16. 414456312 Q= 9.21761357658e-06 dim(res. GG)= 
-17.8764434504 Q= 0.0010604463244 dim(res. GG)= 
-18. 0526932154 0= 1 . 30960846002e -06 dim(res. GG)= 
-21 .2659762563 Q= 6 . 9908080844e-07 dim(res. GG)= 
-21.4084980043 Q= .000232765771 dim(res. GG)= 
-25. 1493689791 Q= 5 . 0551 1356071 e -07 dim(res. GG)= 

Figure 2: Example code demonstrating scripting use of the validation code and its 
output. 
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Furthermore, there are functions to map numerical vectors of length 70 that con- 
tain 35 s + 35 c coefficients of the 'boost' generators of ^(+7) to the language of self- 
dual and anti-sclf-dual four-forms and back. Internally, four-form coefficients are 
represented as Python dictionaries (hash tables). The f ormat_abcd_to_string() 
and f ormat_abcd_to_latex() functions can be used to bring these into human- 
readable form. 

The function abcdJrom_v70() maps a 70- vector to a four- form Python dictio- 
nary. Its inverse is given by v70_f rom_abcd() . Normalization conventions are such 
that if a 70-vector v has a single non-zero entry, then the non-zero entries (192 of 
them) in the corresponding four- form dictionary are all ±1. 

2.3.4 The e7_potential.py module 

This component contains the definitions needed to calculate the supcrgravity scalar 
potential. The most important function is e7_potentialjfrom_v70(), which com- 
putes and returns the potential to numerical accuracy. As this function can also 
return the most important intermediate quantities used in the calculation, such as 
the tensors A\ and A2, when called with an additional dictionary parameter info 
to store data in, the other auxiliary functions defined here - such as e7_A12() - are 
usually not needed. 

Based on the potential function, the function e7_test_gradient() uses finite 
differences to numerically calculate an approximation to the potential's gradient as 
well as its length at a given point. This provides a direct test for stationarity and 
hence allows the numerical validation of solutions. 

2.3.5 The e7_vacua.py module 

This module provides just a single constant, the dictionary vacua_v70 that contains 
70- vectors for all the known stationary points in the scalar potential. Some solutions 
are associated with named keys, such as 'G2', and all of them can be addressed with 
keys of the form '#17', corresponding to the list given in [T3] . 

This module also provides a function e7_verify_vacua() that re-does the au- 
tomatic gradient check for all known stationary points. This function uses the 
function e7_test_gradient from the e7_potential .py module. 

2.3.6 The e7_analysis .py module 

This module provides additional functions to analyze the properties of stationary 
points. The most important ones are: 

• e7_test_stationarityjQ() - to verify whether a point satisfies the Q a bcd + 
^il e abcdefghQ e ^ 9h = stationarity condition. This provides an independent 
test in addition to the direct numerical evaluation of the gradient provided by 
e7_test_gradient () . 

• e7_residual_gauge_group_generators() - to automatically determine a ba- 
sis of the subgroup of spin(8) that is left unbroken by the point in question. 

• e7_susy_A2() - to numerically determine the number of unbroken supersym- 
metries via the eM2 i J = condition. 

• e7_Al_eigenvalues() - to numerically determine the eigenvalues of the A\ 
tensor. In addition to these eigenvalues E, the values — 6\E\ 2 /P are returned 
as well. This gives a second test for unbroken supersymmetry. 

• e7_scalar_massmatrix() - to numerically calculate the re-scaled scalar mass 
matrix given by formula (2.25) in [5|: 

96e- l g- 2 L M = - (|?(V) + ^A 2 m npq A 2m n A X m Z ijkl 



n a mni a j a mij An | ^ yr^kHpq 

K)A 2 k A 2 mnl ~ ~^2n A 2 „ikl ^iivo^ 



i2fc Si 2 " mnl ~ ^^2n "^-2 mkl \ ^ijpq 2 

o 

_i_ A i A jkP-^mnpq-^ 

~r — /i2 mnp A 2ij ^ ^ijkl- 

One normally will want to rc-scalc this matrix by multiplying with — 6/P 
(V being the value of the potential) in order to obtain masses in AdS-units. 
After such re-scaling, the Breitenlohner-Freedman stability bound [M] is — (d— 
l) 2 /4 = -9/4 = -2.25. 

• e7_A12_masses2() - to calculate the fermion masses-squared according to 
formula (2.14) in 0. 

3 Outlook 

Deeper numerical investigations than those presented in [13j suggest that the num- 
ber of yet unpublished stationary points in the scalar potential studied here as well 
as related ones may be fairly high. A number of these solutions may easily turn out 
to be sufficiently interesting to warrant a deeper investigation. One might think 
that the idea of using numerical techniques to find approximate solutions should 
be of very limited use in particular for those solutions that break the gauge group 
(almost) completely: Studying a stationary point numerically in order to get ideas 
that may guide a fully analytic approach seems to be a fairly unrewarding exercise 
unless it has a sufficient degree of residual symmetry. However, this is not neces- 
sarily the case. First of all, it is possible to eliminate almost all of the arbitrariness 
in the raw numbers that is associated with the freedom to apply a spin(8)-rotation 
to a given solution, using the techniques that were employed in [151 113] to minimize 
the number of non-zero coefficient vector entries. There is no fundamental obstacle 
to modifying readily existing techniques to use very high precision arithmetics (with 
hundreds of digits) rather than machine-precision 16-digit arithmetics in order to 
determine parameters very precisely. This then allows one to utilize powerful in- 
teger relation algorithms such as in particular PSLQ [T7| to systematically 'guess' 
even very complicated analytic expressions. To demonstrate this, let us specifically 
consider a number such as the coordinate Z( 2 ) from table A.l in [5J: this is given as: 

z (2) = ^l-^ I V / 2+^)(3 + V3-3 1 /4Vl0) (i) 

To 300 digits, the imaginary part of Z( 2 ) is: 



Imz (2 ) ra -0.209269477042954393209112597384263059035631424 
839824728046049646351961493949022694957040180268 
460854129215451425779247977077536023519697290138 
356018655716825868238074107476584372107488197263 (2) 
113362737435710197764338927716172598148291052464 
018329383941096143123970705213689657726314686822 
409336284550127068 

If we truncate this number to 200 digits and use the PSLQ algorithm to then find 
a vector of integer coefficients for integer powers of this number so that the sum is 
closer to zero than 10 -170 , the PSLQ algorithm produces the following polynomial 
coefficients: 

[-64,0,704,0,-240,0,32,0,-1] (3) 



One could say that this analytic expression allows us to make a testable 'pre- 
diction' for the next 100 digits. In fact, using Newton's algorithm to find the zero 
of interest of this polynomial reproduces Im z/ 2 ) to 100 digits beyond the previous 
200. This is demonstrated in detail the example code file e7_example_mpmath.py. 
In order to run this piece of code, the 'mpmath' Python extension [53] has to be 
installed. 

It seems hence very likely that even for new stationary points that can be ex- 
pected to require even more complicated analytic expressions, obtaining - and most 
likely even rigorously proving - analytic results can be fully mechanized. 

The codebase that has been presented in this work (a) correctly reproduces the 
known properties of the 'classical' non-trivial stationary points, (b) should be simple 
enough to be readily validated against the formulae given in |13j , and (c) should in 
the future help to considerably reduce the effort required to check claims about the 
properties of new stationary points. 
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